LPTh-Ji 10/004 



Electroweak Phase Transition in the U( 1 ) '-MSSM 



Amine Ahriche 1 and Salah Nasri 2 

1 Laboratory of Theoretical Physics, Department of Physics, University of Jijel, PB 98 Ouled 

Aissa, DZ-18000 Jijel, Algeria. 
Department of Physics, UAE University, P.O. Box 17551, Al-Ain, United Arab Emirates. 
Emails: ahriche@univ-jijel.dz, snasri@uaeu.ac.ae. 



2 



Abstract 

In this work, we have investigated the nature of the electroweak phase transition in the U(l) 
extended minimal supersymmetric standard model without introducing any exotic fields. The 
effective potential has been estimated exactly at finite temperature taking into account the whole 
particle spectrum. For reasonable values of the lightest Higgs and neutralino, we found that the 
electroweak phase transition could be strongly first-order due to: (1) the interactions of the 
singlet with the doublets in the effective potential, and (2) the evolution of the wrong vacuum, 
that delays the transition. 

Keywords: baryogenesis, electroweak phase transition, extra singlet, extra gauge boson. 
PACS: 12.60.-i, 12.60. Jv, 14.80.Cp, 12.15.-y. 



1 Introduction 

The matter-antimatter asymmetry in the Universe is observed to be n&/n 7 ~ 10 -10 [lj. If 
this asymmetry is to be explained by microphysics rather than initial conditions, then there 
must be processes occurring in the early Universe that violate baryon number and CP and 
which occur out of thermal equilibrium. It appears that the standard model (SM) satisfies 
all three conditions [21 |3]; the baryon number is not conserved at quantum level due to the 
B + L anomaly [4J, a CP violation source does exist in the quark sector (CKM matrix), and 
a departure from thermal equilibrium could in principle be achieved through a strong first- 
order phase transition [3J. However, detail calculations show that the SM fails to generate the 
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observed baryon asymmetry due to the smallness of the CP violation effect and the weakness of 
the electroweak phase transition (EWPT) [5]. 

In gauge theories, a first-order phase transition takes place if the vacuum of the theory does 
not correspond to the global minimum of the potential. Since it is energetically unfavored, the 
field changes its value to the true vacuum ( i.e., the absolute minimum of the potential). Because 
of the existence of a barrier between the two minima, this mechanism can happen by tunneling 
or thermal fluctuations via bubble nucleation. The electroweak baryogenesis scenario is realized 
when the B and CP violating interactions pass through the bubble wall. These interactions are 
very fast outside the bubbles but suppressed inside. Then a net baryon asymmetry results inside 
the bubbles which are expanding and filling the Universe at the end. 

In the SM, the EWPT is too weak [6] unless the Higgs mass is less than 45 GeV [TJ, which 
is in conflict with present data [8j. But a departure from thermal equilibrium without being 
in conflict with this severe bound on the Higgs mass, is possible when extending the SM with 
additional gauge singlets [9J [TO] , new heavy fermions [TT], an extra Higgs doublet [12], or in some 
supersymmetric extensions of the SM. 

In spite of its success and popularity, the minimal supersymmetric standard model (MSSM) 
with R-parity still has two major problems: the /i-problem [I3J and the fast proton decay due 
to dimension 5 operators [2]. A natural solution to these problems would probably require 
the extension of MSSM by a new mechanism or a new symmetry. The £/(l)'-extended MSSM 
(USSM, UNMSSM or UMSSM) [15] is a straightforward extension of the MSSM with a 
nonanomalous TeV scale Abelian gauge symmetry. This simple enlargement of the gauge sector 
is well motivated in string construction [IB] , in grand unified theories [TTJ such as SO (10) and Eq, 
in models of dynamical symmetry breaking [12] and little Higgs models [TO] . The /i-problem and 
the dangerous dimension 5 operator can be solved naturally with an appropriate U(l') charge 
assignment. Furthermore, these U(l)' models can provide a new candidates for dark matter 
that are not excluded by direct dark matter searches and with interesting signatures at colliders 

[221 EU- 

In the MSSM, the EWPT could be strongly first-order if the light stop is lighter than the 
top quark [22]. In the singlet extended MSSM [23], such as the NMSSM [21], the EWPT get 
stronger easily for a large range of parameters [25]. In gauge extensions of the MSSM, such as 
UMSSM, the EWPT is also strongly first-order but with the price of introducing 3 new extra 
singlet scalars [26], or by adding new extra heavy singlet fermions [27] . 

The main reason that makes it less easier to have a strong EWPT in the UMSSM compared to 
NMSSM , is that the former contains a new gauge interaction which results in strong constraint 
on the mixing between the SM gauge boson Z and the new one Z' [2H] 

2M 2 ZZ ,/(M 2 Z , ZI - M 2 ZZ ) < lO" 3 , (1) 

and the bound on the heavy Z' mass [29J 

M z > > (500 - 800) GeV, (2) 

which implies serious constrains on the vacuum expectation value (vev) of the singlet and the 
new U(l)' gauge coupling g'. 

However, both models have similar form for the scalar potential, where the singlet can play the 
same role during the EWPT dynamics. In this type of mode0, the singlet vev within the wrong 
vacuum could be nonzero, i.e., (S) = i^0, and therefore, is temperature dependant during the 
EWPT dynamics. This feature could delay the EWPT, i.e., lowers the critical temperature, and 
enhances the parameters that define the strong first-order phase transition criterion [50] : 

v(T c )/T c >l, (3) 
1 Similar remark holds for models with a singlet like [51 110). 
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where T c is the critical temperature and v (T) is the temperature dependent scalar vev. 

In this work, we will investigate the possibility of getting a strong first-order phase transition 
within the minimal gauge extension of the MSSM, UMSSM without adding any new field beside 
the usual singlet. 

This paper is organized as follows: in the second section, we give a brief review of the UMSSM 
model, define the effective potential and discuss different constraints on the parameters. After 
that, we discuss the EWPT dynamics and show how to get a first-order phase transition. In the 
fourth section, we discuss our numerical results. Finally, we summarize our results. The different 
field-dependant masses used in the estimation of the effective potential are given in Appendix A. 

2 The UMSSM model 

The If (l)'-MSSM (or UMSSM) is based on the gauge group G = SU(3) C x SU(2) L x U(1) Y x U(l)' 
with the couplings g$, g 2 , g\ and g', respectively [15] and the superpotential is given by 

W = \SeijH\Hi + YueijQWHi + Y D e ij Q i D c H{ + Y L t % jVE c El (4) 

where is the antisymmetry 2x2 tensor, Yjy, Yd and Yl are Yukawa couplings, and A is a 
coupling constant in which A (5) generates the //-term in the MSSM. The particle content of 
this model is given by the left-handed chiral superfields L ~ (1, 2, —1/2, Ql), E c ~ (1, 1, 1, Qe), 
Qi ~ (3,2,1/6,Q Q ), U c ~ (3,1, -2/3, Q v ), D c ~ (3, 1, 1/3, Qd), #1 ~ (1,2,-1/2,^0, #2 ~ 
(1,2, 1/2, Q 2 ) and S ~ (1,1, 0,Qs), where the 1/(1)' charges, Q/s, are model dependent. For 
instance, in a class of Eq gauge models, the group can be broken in two steps to its 5*0(10) and 
51/(5) subgroups: 

E 6 -»■ 50(10) x 1/(1)^ -»■ 51/(5) x U(l) 4 , x U(l) x , 

and in this case, the 1/(1)' generator is given in terms of the U(l)^ and U(l) x generators and 
mixing angle 9e 6 [E] 

Of = cos 9 Efi + Q x sin 9 Ee • (5) 

Although the SU(3)c X SU(2)i x U (l)y x U (1)' model is motivated in the Eq framework, we will 
not single out a particular charge pattern and we only require that the model be anomaly free. 
However, with the above particle content, one can easily show that the invariance of the Yukawa 
terms in the superpotential under U{1)' and the absence of the SU(3)c — SU(3)c — U{\)' anomaly 
implies that Q s = 0, and the // problem arises again. Thus, anomaly cancelation requires exotic 
representations beyond those of the MSSM. The simplest extension for the anomaly to vanish is 
to assume three generations of heavy (few TeV) vectorlike pairs of chiral fields Ki and K? which 
transform as (3, 1) and (3, 1) under SU(3)c X SU(2)l and opposite hypercharges [3Tj . 

2.1 The effective potential 

In case where both squarks or/and sneutrinos do not develop vevs, the scalar potential is a 
combination of the so-called D, F and soft terms, which are given by 

Vb = {H+H 2 - H+H,) 2 + I \H+H 2 \ 2 + ^ + Q 2 H+H 2 + Q s \S\ 2 f , 

V F = \\\ 2 {\t l3 H\Hi\ 2 + |5| 2 [H+H, + H+H 2 ]} } 

= m\HtH x + m 2 2 H+H 2 + m| |5| 2 + {A x Se i:j H\Hi + h.c}. (6) 
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Here mf, m|, m| and A x are usually called the SUSY soft parameters. The charges Q's should 
be chosen in such a way that the anomaly cancelations are ensured [T5J [20] . 

The structure of the tree-level potential (Q seems to allow two explicit relative CP violating 
phases between the doublets and the singlet. However, gauge invariance dictates that the only 
allowed phase in the potential, is that of the combination SeijH[Hl, since gauge rotations can 
be used to set the charged Higgs vev (H^ = 0, and the condition {H 2 ) = implies that 
the physical charged Higgs is nontachyonic (M^ ± > 0) [T5]. Therefore, the general form of the 
ground state could be written as 

/ \ (■■> ( 1 \ (',.( " 



Here v = (v\ + v 2 ) l l 2 = 246 GeV, and 9 is the relative CP violating phase between the singlet 
and the doublets; and 9 its value at the ground state. The tree-level scalar potential in terms 
of the neutral components is given by 

V (h, h 2 , S) = [(gl + gl) /4 + g' 2 Q\] hj/8 + [(gf + g\) /4 + g' 2 Q 2 2 ] hj/8 + g' 2 Q 2 S 4 /8 

+ [A 2 - {gl + gf) /4 + g' 2 Q x Q 2 \ h\h 2 2 /A + [A 2 + g^QxQs] h\S 2 /A 

+ [A 2 + g' 2 Q 2 Q s ] h 2 2 S 2 /4 + m\h\/2 + m 2 2 h 2 2 /2 + m 2 s S 2 /2 + {A x Sh x h 2 e i& jl^ft + h.c}. 

(8) 

The neutral scalar sector in this model contains three CP-even scalars and one CP-odd scalar. 
The tadpole minimization conditions at the ground state along the CP-odd scalar forces the 
relative CP phase 9, to be canceled by the phase of the parameter A\ in which must be 
taken as A\ = ± \ A\\ e~ td ° . This makes the ground state at tree- level independent of this relative 
phase, and the phase S in <^ is just 5 = 9 — 9 . It is clear that the dependence of ([8]) on the 
phases disappears at the ground state. However, when considering the one-loop corrections in 
the tadpole minimization condition along the CP-odd scalar, the tree-level phase cancelation is 
no longer valid at one-loop, and the argument of Ax differs slightly than —9 Q (or —9 Q + ir). 

The structure of the tree-level potential ([6]) implies the presence of a mixing between the 
submatrices (CP-odd scalars, Goldstone bosons, and CP-even scalars) at tree-level for any values 
of the scalar fields except at the ground state, however this mixing is proportional to sin 5. This 
will ensure the CP conservation at tree-level in the scalar sector. Indeed, in addition to the phase 
appearance in the scalar sector due to the deviation of Arg(A\) from — 9q, the CP violation effect 
could be seen in this model at one- loop [32], and also through the dependence on 9 and not 5 
can be seen in the superpartner masses (1231) . fT26|) and fl27|) . 

The parameters mf, m 2 , and m| can be eliminated by taking (v±, v 2 , v x e l9 °) as a minimum 
of the effective potential. The vacuum stability of ([8]) requires the condition: 

A 4 - A 2 [{gl + gl) 12 + 2^' 2 (Q 2 S - Q 1 Q 2 )] + \g< 2 [g\ + g 2 2 ) Q 2 S > 0, (9) 

where the U{1)' gauge invariance condition Q\ + Q 2 + Qs = is taken into account. 

In models that include singlets, like the present one, one needs to be careful about whether 
the minimum (v\, v 2 , v x e ld °) is the absolute one. This can be checked by comparing the effective 
potential value at this minimum with its value at the wrong vacuum, which could be, in singlet 
models, (0, 0, x) rather than the origin (0, 0, 0). The wrong vacuum can be defined as the mini- 
mum of the effective potential in the direction where all the SU(2) doublets vanish. At tree-level, 
it is given by 

M 2 = vl + [Q./Qs + \ 2 /g' 2 Q 2 s ] v\ + [Q 2 /Q s + \ 2 /g' 2 Q 2 s] v 2 2 + -^—v.v,. (10) 
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The value of x in fllOp should be relaxed by including the one-loop corrections given below in Eq. 
(ITT]) . In the case where the right hand side of (fTOl) is positive, the wrong vacuum (0,0, x) does 
exist, otherwise it should be (0, 0, 0). In both cases, (v\, v 2 , v x e °) must be the absolute minimum 
for the effective potential. This condition imposes additional constraints on the parameters of 
the model. 

The one-loop effective potential at zero temperature is given in the DR scheme by [33] 

V 11 (hi, h 2 , S) = V (hi, h 2 , S) + J2^ ( lo S - I) . (11) 

where mi (hi, h 2 , S) are the field-dependant masses, which are given in Appendix A, A is the 
renormalization scale, which is taken to be A = v = 246 GeV, and rii are the fields multiplicities: 
n w = 6, n z = n z > = 3, n h o = n h o = n s = n A o = n G = 1, n t = -12, n- tL = n~ tR = 6, = -2, 

= —4, where A , G, ti,R, X an d C denote the CP-odd Higgs, Goldstone boson, left- and 
right-handed squarks, neutralinos and charginos, respectively. 

The thermal corrections to the effective potential can be computed using the known techniques 
|34j . The one-loop effective potential at finite temperature is given by 

V eff (h, h 2 , S, T) = V 1 - 1 (h u h 2 , S)+T 4 J2 n * J B,F K (h, h 2 , S) IT 2 ) ; (12) 

i 

with 

1 r°° 

J B ,F(a) = — x 2 log(l =Fexp(-v / x 2 + a)), (13) 

rii are given above and mf (hi, h 2 , S) are given in Appendix A. In this work, to take into account 
all the heavy and light fields, we will evaluate this integral numerically. 



2.2 The parameters 

In this model, we have many parameters, some of which are free like: g', A, v x , tan/3 = vi/v 2 , 
and the soft terms: itiq, mu, A t , A\, M 2 , Mi, and M[; and others that are fixed by a measured 
physical quantities such as : gi, g 2 , v and y t ; or can be conditions like the elimination of m\, m\, 
and m| in 

In scanning the parameter space of this model, we take into account the constraint Qi + 
Q2 + Q.S = 0, conditions from the minimization of the potential, the perturbativity of the quartic 
couplings in ([8]), and the vacuum stability ([9]). Another constraint could be derived from the 
upper bound on the mixing between the gauge boson Z and the new one Z' ([T]), and the lower 
bound on the Z' mass (T5]). The condition (TJ]), could be achieved by considering relatively large 
v x or large (g'Q 1 )- The condition ([1]) could be fulfilled if the mixing term M\ z , is vanishing, i.e., 

Q 1 = Q 2 tan 2 0, (14) 

which leads to a fine tuning in the values of Qi )2 and tan/3. The second possibility is making 
M\, z , » M zz , M\ z , , which roughly means 

g'\Qs\ v x > (500 - 800) GeV. (15) 

In our search for the parameter's space that fulfills the strong first-order phase transition 
criterion, v (T c ) jT c > 1, we will focus on the two following regions : 

(1) Moderate values for the parameters Qi ;2 and tanfi, where (pQ) is nearly satisfied. In this 
case, the singlet vev v x can be of order v or relatively smaller. 
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(2) The two terms M\ z and M 2 ZZ , in (|22p . are suppressed with respect to the mass term 
M\, z ,. In this case, the values of U'(l) charge and the vev of the singlet, Q s and v x , must be 
large enough. 

We should also distinguish between the two cases where the minimum (0,0, x ^ 0) of the 
potential does exist or not. The second case could be ensured by choosing the parameter A\ in 
(fTUj) to be extremely negative, then x does not exist. This is easier to satisfy in region (1). In 
both cases, the condition §2§ is, almost, automatically fulfilled within ([I]). 

Also, one needs to know the effect of the phase 9 on the EWPT strength. At tree- level ([H]), 
this phase is not relevant due to the choice of A\, but the thermal corrections depend on this 
phase through all the fields masses expect gauge bosons and top quark. Therefore, we consider 
the values < 9 < n/3. 

The mass parameters: mq, mu, A t , M 2 , Mi, and M[ appear at one- loop level in the effective 
potential, therefore we expect that their role is less important in the EWPT dynamics. Indeed, 
these parameters will change the field masses and therefore could make (vi, v 2 , v x e td °) a local 
minimum instead of the global one through the corrections in ffTTj) . However, the parameters: 
g', A, v x and tan/3, as well the charges Q's that appear multiplied by g', seem to play a very 
important role for the strength of the EWPT. Therefore we focus on these parameters while 
fixing g' = g\, vtlq = mjj = 1 TeV, and choosing different values for A t , M 2 , M 1; and M[. We 
will allow the parameters A, v x , Qi t 2, A\, and 6q to vary randomly within the intervals : 

0.001 < A < 0.5, 0.5 < v x /TeV < 4, 

1 < tan/3 < 20, —1.2 < A\/TeV < 1.2, (16) 
-4 < Q 1>2 < 4, < 9 < tt/3. 



3 The electroweak phase transition 

Due the condition (TjQ) and (j2J), there could exist a hierarchy between the vev of the singlet and 
those of the doublets, i.e., v x » vi^- In this case, the gauge symmetry could be broken in 
two step. However, in the case where the mixing is extremely suppressed as in region (1) (jl4|) . 
the singlet vev could be low as ~ 500 GeV. In this case, the gauge symmetry could be broken 
just in one step. In the two-step symmetry breaking case, one notices that above a certain 
high temperature, the singlet vev was zero v x (T>>) = 0. This can be seen by putting Vi )2 = 
in (llip . and taking only the thermal correction of Z'. Then at lower temperatures, it is not 
sure that the system moves directly from (0,0,0) to (vi, v 2 , v x e t0 °), or via an intermediate step: 

(0,0,0) (0,0, x) -4 (vi,V2,v x e ie °). This will depend, in general, on the theory parameters, 
especially the value of the singlet vev. If it is comparable to the EW vev v, the phase transition 
could occur once, however if it is much larger than v, the phase transition will occur in two steps. 

Since the singlet dynamics does not affect the SU{2) sphaleron processes, we will not be 
interested in distinguishing between the one- and two-steps symmetry breaking. We will treat 
our field dynamics using the effective potential where the singlet is replaced by a temperature- 
dependant vev, i.e., 

V(/ii,/i 2 ,T) = V eff (h u h 2 ,x(T),T). (17) 

At higher temperatures, the effective potential admits only one minimum where (T) = 0. 
As the Universe cools down, the effective potential acquires a new minimum (T) ^ 0; but it 
is not the absolute one. At such temperature, the critical temperature T c , the two minima get 
degenerate 

Veff K, v c 2 , v c x , T c ) = V eff (0, 0, x c , T c ) . (18) 

Below this temperature, the new minimum becomes the absolute one, and the system has to 
move from the old (false) vacuum to the new (true) one. In the case where a barrier does exist 
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between the two minima, this transition has to occur via tunneling trough bubbles nucleation 
at certain points, which expand and the whole space by the new vacuum (T) 7^ 0, i.e., the 
symmetry is broken. 

If the effective potential ()12p is expanded as powers of m/T (in the limit m « T) in a 
similar way to the SM in the doublets direction, the two leading terms are of order h 2 T 2 and h 3 T. 
The first term determines the temperature when the barrier between the potential two minima 
disappears, while the second term is relevant to the strength of the first-order phase transition. It 
turns out that a resummation of the so-called Daisy diagrams [35J, leads to a contribution of order 
h 3 T which exactly cancels the contribution from certain particles (e.g. Higgs and longitudinal 
gauge bosons). Therefore this screening effect could weaken the strength of the EWPT. For that, 
we will check the importance of this effect by using a resumed effective potential where (Tl2|) is 
modified by replacing the bosonic field-dependant masses by their thermally corrected values 
[55] . The thermal corrections to the bosonic masses are given in Appendix B. 

The B + L anomalous interactions [I] that violate the baryon number do not have the same 
rate in the symmetric and broken phases ( i.e., at both sides of the bubble wall). In the symmetric 
phase, this rate behaves like ~ T 4 [57], and is suppressed as exp(—E Sp /T) [38], in the broken 
phase, where E$ p is the system static energy within such field configuration called the sphaleron 
[59] . Therefore any generated baryon number at the symmetric phase will get erased , unless these 
interactions are switched off in the broken phase, which translates to the condition ([3]). In reality, 
the B + L anomalous interactions should be switched off at a temperature Tq < T n < T c , where 
T n is the temperature at which bubbles start to nucleate, and To is the temperature at which 
the barrier between the two minima completely disappears. Then the condition in ([5]) should 
be fulfilled at T n , but, in general, the two values are significantly close, and if v(T c )/T c > 1, it 
is necessarily satisfied at T n . Indeed, in some earlier works (e.g. [I0J E]), the condition (J5J) at 
To was used in defining the phase transition. This means that the anomalous B + L violating 
processes should be switched off at the end of the phase transition, i.e., the Universe is filled by 
the true vacuum by the expanded bubbles. The value T , usually called the lower metastability 
temperature or lower spinodial decomposition point [UJ, can be determined using the Jacobian 
det[d 2 V(h, T )/dhidh 

i]h=o — 0- However, in our work, it is safer to consider the condition ([5]) 
at T c , that is defined in (120]) . 

Since the singlet field does not play an important role in the sphaleron processes [TO] . the 
criterion of a strong first-order phase transition in our case is given by 



In the general case where the relative phase 6 7^ 0, the field ground state at nonzero tem- 
perature should be written as {vi} i=14 = (vi, v 2 , v x cos9, u^sin^) instead of (vi,v 2 ,v x ) and the 
relative phase. These 4 variables should be treated independently when looking for vi^, x and 8 
at any temperature T. Then the phase transition could be defined through the equations 



where are the real and imaginary parts of the T-dependent x given in (flOl) . 

4 Numerical results 

In the following figures, we show the dependence of the quantity v (T c ) /T c in (1T9]) on the lightest 
Higgs mass (in Fig. [JJ left) , and on the lightest neutralino mass (in Fig. [JJ right), for a random 
choice of about 10 6 cases in both regions (1) and (2). We find that only about 8% of the 




(19) 



d_ 

dvi 



V eff (Vi,T c ) = 0, V ef f {v U V 2 ,V Z ,V^T c ) = V eff (0,0,X 1 ,X 2 ,T c ) , 



(20) 
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benchmarks fulfill the required conditions, and only 2.5% of the survived benchmarks give a 
strong first-order phase transition. 




Figure 1: The dependence of the quantity v c /T c on the lightest Higgs mass, on the lightest neu- 
tralino mass, and the corresponding critical temperature, taking into account different values of 
the parameters in j[Tb)) is shown. The red points refer to benchmarks from region (1), and the 
green ones to benchmarks from region (2). 

As it is clear from Fig. HJ the EWPT could be strongly first-order in the two regions (1) 
and (2), for a large range of the of the values of the lightest Higgs mass, and lightest neutralino 
mass. These masses are estimated at tree-level, which are expected to be boosted to larger 
values when considering one-loop corrections. One remarks that the phase transition strength 
does depend on the lightest Higgs and neutralino masses. One also remarks that the majority 
of the benchmarks that give a weak phase transition have the value of v (T c ) /T c ~ 0.14 — 0.18, 
and their corresponding critical temperature is between T c ~ 300 — 600 GeV. While all the 
corresponding values of T c to the strong phase transition are T c < 180 GeV. 

The mixing with the singlet will modify the doublets couplings to the fermions, and since 
the lightest Higgs contains a singlet amount, its mass could be (much) smaller than SM bound 
~114 GeV. One remarks from Fig. [1] (right), that the EWPT could be strongly first-order for 
small values of the lightest neutralino masses; m Xl < 15 GeV. Such a light mass of neutralino 
with a spin-independent cross section of order 10~ 5 pb could be a possible interpretation of the 
recent observations of CoGeNT and DAMA/LIBRA |42J. In particular, the lightest neutralino 
can have a non negligible component of the superpartner of the U(l)' gauge boson [15] . 
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The critical temperature is, in general, higher when comparing with minimal SM (~100 GeV), 
the generic value is larger than 300 — 600 GeV. This is a consequence of the interaction of the 
doublets with the singlet that has, in general, a very large vev. However for the benchmarks, 
giving a strong first-order EWPT, T c is relatively smaller than the generic value. In fact, for 
points (1) as T c is smaller than about 180 GeV, the stronger the EWPT becomes. 

In order to understand this point, we take a benchmark from Fig. [TJ and study the dependence 
of the scalar vevs on the temperature T below and just above the critical temperature. We will 
also check how could this behavior be changed with respect to the charges Q's, and other 
parameters like A t , and M[, that appear in the effective potential at one- loop level. Therefore 
we consider the benchmarks in Table-[H 





(a) 


(b) 


(c) 


(d) 


A 


0.0235 


0.0235 


0.0235 


0.0235 


tanj3 


2.0566 


2.0566 


2.0566 


2.0566 


A x 


-208.9569 


-208.9569 


-208.9569 


-208.9569 


Vx 


1173.3560 


1173.3560 


1173.3560 


1173.3560 


Qi 


0.3575 


1.75 


0.3575 


0.3575 


Q 2 


1.6147 


1.5 


1.6147 


1.6147 


0o 


0.5888 


0.5888 


0.5888 


0.5888 


Mi 


857.9520 


857.9520 


857.9520 


1500 


M 2 


398.7435 


398.7435 


398.7435 


398.7435 


M[ 


881.3455 


881.3455 


881.3455 


881.3455 


A t 


274.4727 


274.4727 


1400 


274.4727 


T c 


103.6387 


326.2436 


302.9225 


100.8061 


v c /T c 


2.5836 


0.1940 


0.1976 


2.6356 


m hl 


172.1886 


163.3732 


173.0714 


172.1589 


m Xl 


5.6907 


5.6569 


5.6907 


6.2573 



Table 1: The values of the parameters used to study the scalar vevs dependence with respect to 
the temperature are shown. The mass- dimension parameters are given in Ge V. 

From this table, it is clear that the EWPT strength is extremely sensitive to the charges Q's 
which represent in a way the strength of the U(l)' interactions since the Q' charges always appear 
multiplied by g'. However, the soft parameter A t and Ml that appear in the effective potential 
at one loop could have an important effect on the strength of the EWPT for some particular 
choices of the parameter space. The dependence of the EWPT strength on the relative phase 
#o, could be seen by taking changing its value of the benchmark (a) in Table-HJ We find that its 
effect is extremely negligible [difference in the value of v (T c ) /T c is < 0.1%]. 

In Fig. 121 we show the dependence of the ground state on the temperature below the critical 
temperature for the benchmark (a) in Table-HJ and its modifications (b), (c) and (d). One 
remarks that the common feature between all these cases is that the dependence of the singlet 
vev on the temperature is very weak around and below the critical temperature. This could be 
understood due to the fact that the singlet vev is much larger than this temperature [x(T « 
x) ~ 3^(0)], and also, due to a possible two-stages phase transition realization, where the first 
stage (0,0,0) — > (0,0, x(T)), takes place around the temperature T* ~ x; therefore at lower 
temperatures T << T*, the singlet vev becomes almost temperature independent. 

Another important remark, is that the Higgs vevs for this benchmark [Fig. [2]-a] are increasing 
when the Universe gets cooled unlike the SM [6J, or MSSM [41], and similar remark for (d). In 
cases (b) and (c), the Higgs vevs decrease, but with slower pace than in the SM or MSSM, which 



9 




Figure 2: The dependence of the scalar vevs on the temperature T is shown. These quantities 
are scaled by their zero temperature values. The solid curves refer to the broken phase, where 
(v\,V2,v x ) represents the global minimum. While the dashed ones refer to the symmetric phase, 
where the global minimum is (0,0,x(T)). 

makes the critical temperature larger than in the SM. The modified cases (b) and (c) belong to 
the majority of the benchmarks mentioned before. 

The effect of the increasing doublets' vev values at high temperatures with respect to their 
zero temperature values has been mentioned in a similar work [27]. This behavior and the slow 
decrease of the doublets vevs with respect to the temperature, is a consequence of the interaction 
of the singlet with the doublet Higgs. These interactions have the effect of relaxing the shape of 
the potential in the direction of the doublets, and therefore enhancing the ratio in ( II 9p . This is 
a common feature for models with singlets [TO] . 

From ()20p . the phase transition is defined as the degeneracy of the two vacua. Unlike the 
SM and MSSM, the wrong vacuum (0, 0, x(T)), is evolving with respect to the temperature, and 
could be reached through a phase transition at very high temperature. Therefore, its evolution 
could be a very important factor which can strengthen the EWPT. It is defined from the effective 
potential as the local minimum 

V r (T) = K//(0,0, x(T),T) , A K//(0 ,0,x(T),T) = 0, (21) 

which was a global one before the phase transition, i.e., above the critical temperature. We need 
to check that the evolution of (12"T|) is responsible of lowering the critical temperature in (a) and 
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(d), which does not play the same role in (b) and (c). In Fig. [31 we show how does the evolution 
of the wrong vacuum with respect to the temperature decrease the critical temperature for the 
benchmarks (a) versus its modification (b). 




Figure 3: The evolution of the effective potential values at the minima (0,0, x) and (v\,v-i,v x ) 
with respect to the temperature T in units of v A is shown. In this temperature range, the value 
of the effective potential at (0,0,0) is about 5.879^6.620 for benchmark (a), and 7.322^8.131 
for benchmark (b) in units of v A . 

From Fig. [31 one remarks that the minimum (vi,V2,v x ) becomes the true one at low tem- 
perature for case (a), and at large temperature for case (b) due to the evolution of the wrong 
minimum (0, 0,x). This can be identified from the intersection of the thermal effective poten- 
tial at (vi,V2,v x ) and its value at (0, 0,x). If one completely ignores the minimum (0,0, x), the 
EWPT, which is defined in this case by the intersection of this value with effective potential 
at (vi,V2,v x ) and its value at (0,0,0), will take place at very high temperatures, normally of 
the order ~ 0(v x ). It is clear that the wrong vacuum evolution has the main role in lowering 
the critical temperature in case (a) compared with case (b). It might seem that the EWPT is 
less sensitive to some parameters such as ttlq, my, A t , M2, Mi, and M{, since they appear in 
the effective potential at one-loop. However, we have shown in Table-ID and in Fig. [3] that the 
EWPT dynamics is sensitive to these inputs more than expected. In fact, if one adds a small 
perturbation to these parameters, then the wrong vacuum, (0, 0, x), becomes the absolute one at 
zero temperature, through the one- loop corrections in (fiT|) . which will rule out our benchmark. 

In order to estimate the screening effect on the EWPT strength, we consider the modified 
effective potential [36] by replacing the longitudinal gauge bosons and squarks masses in (TT2|) 
by their thermally corrected expression (1301) and (1311) in Appendix B. We ignore the scalar 
contributions because they are less relevant to the EWPT dynamics. One can distinguish two 
different types of the benchmarks in Fig. [U the first variant: the benchmarks with increasing 
doublet vevs with respect to the temperature [like the benchmark (a) in Fig. [2], and the second 
variant are the benchmarks with decaying doublet vevs with respect to the temperature like most 
of the benchmarks in Fig. [TJ 

For the second variant, it is expected to have a similar behavior as in the (MS)SM-like models 
where the screening effect weakens significantly the strength of the EWPT. The benchmarks 
whose generic values of v c /T c in the interval 0.14 — 0.18 get reduced to 0.10 — 0.12, that is about 
30% effect. In fact, even the benchmarks that give a strong first-order EWPT but with (slowly) 
decaying doublet vevs with respect to the temperature become weak EWPT benchmarks after 
including the daisy- diagrams contribution. 
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Without including the daisy-diagrams, the first variant correspond to benchmarks giving a 
strong first-order phase transition whose critical temperatures are less than T c < 180 ~ 200 GeV. 
Once we include the daisy diagrams contribution, we find that the doublets vevs become sharply 
increasing with respect to the temperature, and this results in a very low critical temperature 
(T c ~ 32 — 36 GeV), and very strong (v c /T c ~ 7 — 9) EWPT. This unusual effect is due to the 
breakdown of the approximation in which one neglects the mT term in the expression of the 
thermal mass. In other words, one should expect that at low temperatures the thermal mass 
is approximately given by couplings x zero temperature mass instead of the T 2 term [35J. For 
perturbative couplings, this is smaller than the zero temperature mass of the particle and can 
be neglected. Thus, unlike the second variant, we expect that the screening effects for the first 
variant to be negligible and do not significantly reduce the strength of the phase transition found 
earlier (which is first-order). To estimate how small the screening effect is, requires taking into 
account the exact thermal corrections in the effective potential which is beyond the scope of this 
paper. 

It is worth noticing that large values of v (T c ) /T c [i.e., larger than v (T c ) /T c > 3.5] can be 
easily obtained. This corresponds to a severe suppression of the sphaleron (B + L violating) 
processes inside the bubbles, and not necessary to the freezing of the Universe in the wrong 
vacuum. The decay of the wrong vacuum is related to the bubbles' dynamics, not to the (non- 
efficiency of the B + L violating processes. This point requires a special careful investigation to 
put constraints on the theory parameters from the fact that the wrong vacuum must decay into 
the true one. 

5 Conclusion 

In this work, the nature of the electroweak phase transition within the minimal U(l) extension of 
the MSSM (UMSSM) without including exotic particles, has been investigated. We found that 
the EWPT could be strongly first-order for a large range of the lightest Higgs and neutralino 
masses, without adding extra singlet scalars or fermions. We evaluated the effective potential 
at one-loop taking into account the whole particle spectrum, and its temperature-dependant 
corrections were estimated exactly using the known techniques. 

We found that the strength of the EWPT could be enhanced due to two factors: first, the 
interactions of the singlet scalars with the doublets which relax the shape of the effective potential 
in the doublets directions, and therefore, lead to a large value for the ratio v (T) /T at the critical 
temperature. The second factor is that the temperature-dependant local minimum, (0, 0, x (T)), 
could play an important role during the EWPT dynamics. It can delay the phase transition until 
relatively low temperatures (even below 100 GeV), which favor the ratio v (T c ) /T c to be large 
enough, without conflicting the usual severe experimental constraints of the SM and MSSM. 

During this dynamics, the doublets vevs could be decaying with respect to the temperature 
but slower than in the case of the SM or MSSM. Another unusual behavior is that the doublets 
vevs could be increasing with respect to the temperature, which leads to a very strong first-order 
EWPT. We found that the inclusion of the ring contribution does weaken the EWPT if the 
doublets vevs are decaying with respect to the temperature (the second variant), and do not 
change significantly the strength of phase transition in the opposite case. 

We also mention that the strength of the EWPT, as well as the reliability of the theory 
benchmarks, are very sensitive to the input parameters that appear in the effective potential at 
one- loop. 
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A Field-dependent masses 



Here, we will present field-dependant masses with the existence of CP phases; and by vanishing 
these phases, we get the CP concerning case. 

Gauge bosons. In this model, the gauge bosons masses do not depend on CP phases. They 
have similar masses as in the MSSM, but the Z boson is mixed with the new Z' 

M 2 W = \gl (hl+hl) , 

M 2 ( i (gf + g j) (hl+hl) \g'\fgT+g\ [Qih\ - Q 2 h 2 2 ) \ 
z ~ z ' Kb'VgT+gliQihl-Qihl) g' 2 [Q\h\ + Q^+Q|S 2 ) J' [ZZ) 

Tops and Stops. The stop masses are the same as in the MSSM with replacing \x but its 
effective value n = \S/ a/2: 



2 1 2 1,2 

m t = 2% h 2 



m 



l + ml + (gi-lgj)(hl-hl)/8 A t h 2 /V2 + \hSe»/2 \ 



* V Ah/V^+XhSe-^^ m l + m 2 t +gj(hl-hl)/Q J' 

Scalars. The spectrum of physical Higgses after symmetry breaking consists of three neutral 
CP even scalars, a mixture between one CP odd pseudoscalar (A ), and two Goldstone bosons 
that are absorbed by Z and Z' respectively, and a mixture between a charged Higgs and the 
Goldstone bosons that are absorbed by W ± . In the ground state hi = (hi), there is no mix- 
ture between the real and imaginary parts of the complex fields of the scalar sector, and the 
squared-mass matrix should be represented in two independent 3x3 matrices in the two basis 
[V2Re(Hf), V2Re(H°), ^Re(S)} and [ v / 2/m(/J 1 °), v / 27m(//°), V2Im(S)]. However, in general 
case where hi ^ (hi), such a mixing does exist and it is proportional to sin S, as shown in Sec. 2.1, 
and this mixing vanishes at the ground state. Then the elements of the scalar field-dependent 
mass-squared matrix in the basis [V2Re(H?), V2Re(H$), y^2Re(S), V2Im(H%), y^Im(H^), 
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y/2Im(S)], are given by 

Mu = (g\ + gl) {3hl - h 2 2 ) /8 + g' 2 Qi (ZQthl + Q 2 h\ + Q S S 2 ) /2 + A 2 (h\ + S 2 ) /2 + m 2 , 
M 2 12 = (A 2 - (gi + gl) /4 + g l2 Q x Q 2 ) h x h 2 + Re(Q)S, 

M 2 3 = (A 2 + g' 2 QiQ s ) hS + Re(Q)h 2 , M 2 4 = 0, M 2 5 = Im(Q)S, M 2 6 = Im(Q)h 2 , 
M22 = (gl + fe 2 ) (3/^2 - hi) /8 + </ 2 Q 2 (3Q 2 /i2 + Qi^? + Qs^ 2 ) /2 + A 2 (/i 2 + S 2 ) /2 + m 2 , 
^23 = (A 2 + 2 ,2 Q 2 Qs) h 2 5 + Re{Q)h u M 2 4 = Im(Q)S, M 2 5 = 0, M 2 6 = Im{Q)h u 
Mi 3 = g' 2 Q s {SQsS 2 + Qthl + Q 2 h 2 2 ) /2 + A 2 (/i 2 + /i 2 ) /2 + m|, M 3 2 4 = Im(Q)h 2 , 
M 2 5 = Im(Q)h u Ml, = 0, 

M 2 4 = (# 2 + gi) (hi - hj) /8 + g' 2 Q 1 (Q 1 h 2 1 + Q 2 h 2 2 + Q S S 2 ) /2 + A 2 (h 2 2 + S 2 ) /2 + m 2 , 
M 2 5 = -Re(Q)S, M 2 m = -Re{Q)h 2 , 

Ml, = (gl + gl) (h 2 2 - hi) /8 + g' 2 Q 2 (Q 2 h 2 2 + Q 1 h{ + Q S S 2 ) /2 + A 2 (h\ + S 2 ) /2 + m 2 , 
M 5 2 6 = -Re(Q)h u Ml, = g' 2 Q s (Q S S 2 + Q x h 2 x + Q 2 /i 2 ) /2 + A 2 (h 2 x + /i 2 ) /2 + m|, (24) 

with Q = A x e< 9 - e °yV2. The masses of the charged scalars are given in the basis (H± , H 2 ) 

M 2 U = ((gl + g 2 2 ) hi - (gl - g 2 2 ) h 2 2 ) /8 + g' 2 Q x \Q x h\ + Q 2 h 2 2 + Q S S 2 } /2 + A 2 S 2 /2 + m 2 , 
M 2 2 = (g 2 2 - 2X 2 )hh 2 /A - Q*, M 2 = (M 2 2 )* , 

M 2 2 2 = ((0? + g 2 2 ) h 2 2 - (gl - g 2 2 ) h 2 ,) /8 + g' 2 Q 2 [Q x h\ + Q 2 h 2 2 + Q S S 2 } /2 + A 2 S 2 /2 + m 2 . (25) 

Charginos and Neutralinos. The chargino masses are the same as in the MSSM with 
replacing /x by its effective value /i e /j = XS. The two chargino X12 masses are given by the 
MSSM formula 



m 2 ± 

Xl,2 



A 2 S 2 /2 + M 2 2 + g 2 2 (hi + h 2 2 ) T { (X 2 S 2 /2 - M 2 + g 2 2 (h 2 , - h 2 2 )) 



+2gl (X 2 h 2 2 S 2 /2 + Mlhl) + 2^f2g 2 2 XM 2 h l h 2 S cos 9 



(26) 



where M 2 is the SU{2) gaugino mass. 

In the neutralino sector, there is an extra U(l)' zino and the higgsino S as well as the four 
MSSM neutralinos. The 6x6 mass matrix reads, in the basis (_£?', B, W3, H®, H 2 , S) 



M x° 



( M[ 




g'Qihi 
g'Q2h 2 
\ g'QsSe~ 





Mi 


-\gihi 
\g1h2 






M 2 
\g2h1 
-\g2h2 



g'Qihi 
-\g\h\ 
\g2h1 


~T2 h2 



g'Q2h 2 
\g\h 2 

-\g2h2 



A h 



g'QsSe* 9 \ 




-f, h2 
J 



(27) 



where Mi and M[ are the gaugino masses associated with U(l) and U(l)', respectively. 



B Thermal corrections to the bosonic masses 

All the bosonic fields acquire thermal mass corrections from the three typical diagrams in (jH-a)-(j3]- 
c). In Fig. @]-c, we show the fermionic contributions to the scalar mass-squared matrix elements, 



14 



while scalar contributions could be evaluated from diagrams in Fig. flU-a) and flU-b). The gauge 
contributions could be deduced from Fig. flU-a) and (jH-b) just by replacing the internal scalar 
legs by gauge ones. The thermal corrections to gauge field masses could be obtained similarly to 
the scalar ones by replacing the scalar external legs by the gauge ones. 



' 1 4 

i i 1 

\ / \ 

_ ^ .» - ^ s ' 




(a) (b) (c) 



Figure 4: The one-loop diagrams that contribute to the scalar thermal corrections. 

When comparing the leading terms from each diagram, T 2 /12, —mT/8ir and — T 2 /24 respec- 
tively, one easily finds that the integral contribution of (jH-b) is less important with respect to the 
others, and therefore one could ignore it when taking the high temperature expansion. The high 
temperature expansion is assumed to be a good approximation for temperatures m/T < 2.2 for 
bosonic and m/T < 1.6 for fermionic masses with an error less than 5%. 

In our model, we are interested in temperatures around the EW scale where the singlet 
has already developed its vev, then all the masses and vertices should be evaluated under the 
condition 

(/nA,S) = (0,0,^e* e °). (28) 

Under this condition, many fields will decouple and therefore do not contribute to the thermal 
corrections like \/2Re(S), \/2Im(S), Z', the stops, the neutralinos and the charginos, and our 
case meets exactly the MSSM case [40]. The internal lines inside the loops should be mass 
eigenstates, then the vertices required to evaluate (jl]) need to be modified by the induced mixing 
between the interaction states and the mass eigenstates under the condition fl28|) . In many mod- 
els, scalar contributions to the effective potential are generally neglected due their less relevance 
with respect to the gauge contributions, therefore, we will not consider their thermal corrections 
here. In what follows, we give the thermal corrections to the scalar, gauge bosons, and squarks. 
Gauge bosons. The thermal correction to the charged gauge bosons W , is given by 

^glT 2 , (29) 
while the corrections to the mass-squared matrix elements in the basis {7^, Z^, B'^} are given by 

Ml = ^g 2 2 (17g 2 + I5gl)/(g! + g\)T\ 
T 2 I T 2 

M 2 l2 = —glg 2 /^gl + g 2 2 + y9M9 2 2 - 9 2 )/{gl + gi), 
Ml = -^gig2g'(2Qi - 2Q 2 - 3Qt)/y/tf + &, 

M 2 2 = ^rgigl/^gl + gl + ^(H^ + 8g 2 2 g 2 + 9g 4 2 )/(g 2 1 + g 2 2 ), 

M 2 2Z = ^Qtgig' + ^g'((Qi - Q2) {Sgf + g 2 2 ) + 3Q t gt)/yJgl + gl, 

Mi 3 = ^9\Q\ + Ql + QtQr), (30) 

where Q^t is the right-, left-handed top charge under U(l)'. One can obtain the thermal cor- 
rections for the MSMM by putting g' = in (1301) . however the result is slightly different than in 
because here we consider the scalar contributions. 
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Squarks. The left- and right-handed squarks are not mixed under the condition ( I28p . and 
their thermal corrections to the diagonal elements are 

m} L = ~/ s T> + \glT* + ^T 2 + \y\T\ 

K = \glT> + + \y\T\ (31) 

where g s and y t are the strong and Yukawa couplings. 
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